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Abstract: 

The Jahn- Teller model with E®f3 electron-phonon coupling and local (Hubbard-like) Coulomb interaction 
is considered to describe a lattice system with two orbitals per site at half filling. Starting from a state 
with one electron per site, we follow the tunneling of the electrons and the associated creation of an 
arbitrary number of phonons due to electron-phonon interaction. For this purpose we apply a recursive 
method which allows us to organize systematically the number of pairs of empty/doubly occupied sites 
and to include infinitely many phonons which are induced by electronic tunneling. In lowest order of 
the recursion (i.e. for all processes with only one pairs of empty/doubly occupied sites) we obtain an 
effective anisotropic pseudospin 1 /2 Heisenberg Hamiltonian Hef / as a description of the orbital degrees 
of freedom. The pseudospin coupling depends on the physical parameters and the energy. This implies 
that the resulting resolvent (z — Heff{z))~^ has an infinite number of poles, even for a single site. Heff 
is subject to a crossover from an isotropic Heisenberg model (weak electron-phonon coupling) to an Ising 
model (strong electron-phonon coupling). 

1 Introduction 

It has been known for a long time that orbital degrees of freedom in a systems with Jahn- Teller cou- 
pling can be described by an effective pseudospin Hamiltonian 12 El ■ The coupling parameter of the 
pseudospin interaction is t^/U, where t is the orbital hopping rate and U is the strength of the on-site 
Coulomb interaction. This description is of great interest because it provides a model to study orbital 
ordering and the possibility of orbital liquids in terms of conventional spin theories. 

In this paper the influence of the electron-phonon coupling strength on the effective pseudospin Hamil- 
tonian will studied in detail. In order to keep the calculations simple only the case of a system with E^(3 
Jahn- Teller coupling is considered, and the electron spin is neglected. A recursive projection formalism 
^ is applied to derive the effective pseudospin Hamiltonian. This approach provides pseudospin coupling 
parameters that depend on the electron-phonon coupling strength. 

The paper is organized as follows. In Sect. 2 the model is dcflned. As a physical quantity the 
resolvent, related to the electron-phonon Hamiltonian, is considered. Its relation with physical quantities 
is discussed in Sect. 2.1. The recursive projection method is briefly described in Sect. 3 and the effective 
pseudospin Hamiltonian, obtained from this method, is presented in Sect. 4. Finally, the crossover from 
weak to strong electron-phonon coupling is studied in Sect. 4.1. 



2 The [3 Jahn-Teller Model 

The Jahn-Teller model describes fermions with pseudospin a =|, J,, coupled to phonons. It is defined by 
the Hamiltonian H = Ht + Hq, where Ht is the hopping term of the fermions between nearest-neighbor 
sites j and j' 

Ht^~tY^ c\^cy„ + h.c. 

<jj'><T=T,i 
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and Ho is a local (Hubbard-like) interaction and a phonon term: 

j 

for dispersionless phonons with energy lvq- For a given ensemble of fermions, represented by integer 
numbers nj^ = 0, 1, the Hamiltonian Hq can be diagonalized with product states 

]J|A/j,nj|,nj4), 
j 

where A/j = 0,1,... is the number of phonons at site j. The state ^^jt) '^ji) is an eigenstate of the 
phonon- number operator bjbj for nji = njj^ = and for nji = nj^ = 1: 

|7Vj, 0, 0) = ATj \Ni, 0, 0) 1, 1) = TVj |iVj, 1, 1). 

The corresponding eigenstates with a single fermion at j are obtained from | A'j, 0, 0) as 



|iVj, 1, 0) = c]^ exp[— ^(6] - 6j)]|7Vj, 0, 0), |Afj, 0, 1) = c]^ exp[^(6] - 6j)]|iVj, 0, 0) (1) 

and 



|7Vj,l,l)=c]^c]jiVj,0,0>. 
Hj l-^j) '^jT) '^ji) is an eigenstate of Hq with energies 

Wo 
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The groundstate of J?o has no phonons. 
2.1 The Resolvent 

In the following the resolvent 



shall be studied. It is directly related to a number of physical quantities. One is linked with the ther- 
modynamic properties of a statistical ensemble governed by the Hamiltonian H through the Boltzmann 
weight at inverse temperature /?. It reads 

= [ {iz + H)-^e'f'^^. (3) 
Jt 27r 

r is a closed contour that encloses all eigenvalues of H . Another connection is with the dynamics of 
quantum states in a system characterized by the Hamiltonian H: The evolution of a state |\E'o) at time 
to a later time t is given by 

|M/t)=e*^*|fo). 

A Laplace transformation for positive time gives with Imz < the resolvent that acts on the initial state: 

POO POO 

/ e-''^\^t)dt= e-"*e'"Ut\^Q) = {z-H)-^\^o). (4) 
Jo Jo 

The return probability to the initial state is obtained from the inverse Laplace transform and reads 

(*o|*t) = :^/ e"'{^o\{z-H)-^\^o)dz. (5) 
27r J_oo 



2 



Using the spectral representation of H with eigenvalues Ej , the expectation value in the integrand reads 

(vI.o|(z-i7)-|*o) = E^f^- 

j 

Inserting this into Eq. |(SJ allows us to apply Cauchy's Theorem to perform the integration. 

In order to evaluate the resolvent a standard procedure is to expanded it in terms of Ht in a Neumann 
series 

(z - H)-' = {z-Ho- Ht)-' = {z~ Ho)-' Y.^Ht{z - Ho)-'Y 

;>o 

and to truncate this series. The poles of any finite truncation are the eigenvalues of the unperturbed 
Hamiltonian Hq. This may be insufficient for a good approximation of (z — H)-'. In the next section a 
recursive approach is applied that avoids this restriction. 



3 Projection Formalism and Continued-Fraction Representation 

Considering a half-filled system, the projected resolvent Po{z — H)-'Pq must be evaluated. It is assumed 
that Pq projects the states of the entire Hilbert TC space to the subspace Hq- The projected resolvent 
satisfies the identity 



Po(z-H)-'Po 



Poiz - H)Po - PoHPiiz - H)-'PiHPo 



(6) 



where Pi — 1 — Pq projects onto the Hilbert space Hi that is complementary to Hq. If H obeys the 
relations 

PoHPi=PoHP2, PiHPo=P2HPo (P2^Pl), 

Eq. © can also be written as 



Po{z - H)-'Po = \Po{z - H)Po - PoHP2iz - H)-'P2HPo 



(7) 



The identity used in Eq. ^ can be applied again to P2{z — H)i'P2 on the right-hand side. This 
procedure can be applied iteratively. It creates a hierarchy of projectors Pk onto Hilbert spaces Hk- It is 
based on the fact that the projector P2j+i is created from P2j-i and P2j as 



H2j+1 — H2j-l\H2j C H2j~l, 



and -P27+2 comes from the relation 



P2j+iH P2j — P2j+2H P2j = Hj^ij and P2jHP2j+i — P2jHP2j+2 = Hjj^i. 



(8) 



In terms of the projected resolvents this construction implies a recursion relation. Using G2j = P2j{z — 
H)2j-iP2j and Hjj+i = P2jHP2j+2 this reads 



G2j - 



z - P2,HP2j - Hj ,,iG 



j,j + l'^2j+2Hj+ij 



2j 



(9) 



It is useful that the Hamiltonian H in Sect. 2 is the sum of two Hamiltonians H = Hq + Hi and that 
one can choose projections such that the following holds: 

(1) Hq stays inside the projected Hilbert space: HqPq = PqHqPq and Po^^^o = PoHqPq. 

(2) Hi maps from H.2j to H.2j+2- 

Hi : H2j H.2j+2, 

where H2j is orthogonal to H2j+2- In the next section these properties will be used to construct an 
effective Hamiltonian. 
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4 The Effective Hamiltonian 



The projection formalism is now applied to the Hamiltonian H = Hq + Ht of Sect. 2. The Hilbert 
space separates into subspaces with a fixed number of fermions with pseudospin f and a fixed number 
with pseudospin J, because H cannot change it with a diagonal pseudospin term. The case is considered 
here where Pq projects onto singly occupied sites with no phonons. Ht is off-diagonal with respect to the 
phonons and changes the number of pairs of empty/doubly-occupied sites (PEDS) by one. Therefore, P2j 
(j ^ 1) projects onto states with j PEDS and any number of phonons. According to this construction, 
the matrix elements of P2jHP2j+2 P2jHtP2j+2 and P2j+2HP2j = P2j+2HtP2j are all non-zero with 
respect to different phonon numbers. P2jP[{)P2j is diagonal in the basis (^3) with matrix elements: 

^oY^Nj- 2{M - j)g^luo + Uj, 
j 

where 2{M — j) counts the number of singly-occupied sites and j the number of doubly-occupied sites on 
a lattice with 2M sites. Thus the recursion relation of Eq. ([TJ becomes 



G 



2j 



T 



2j 



(10) 



with Hj j+i defined in Eq. The recursion terminates on a finite lattice (M < oo) if j — M, since 
at most M PEDS can be created and, therefore, P2M+1 is a projection onto the empty space. G2M is 
diagonal with matrix elements 

1 

z-u;oEiNi-UM- 

This can serve as a starting point for the iterative approximation of Gq. On the other hand, we can 
express Go by G2 through Eq. H10|l and approximate G2 by the diagonal matrix 



G2^[z-iUoJ2 ^Ih + 2(M - l)5V^o - U] 



-1 
2 • 



(11) 



This approximation corresponds with a truncation of all scattering processes with more than one PEDS. 
It is valid for weak hopping, i.e. for t/ujo small in comparison with g/ujQ and U/luq. Subsequently, it will 
turn out that an effective expansion parameter for the truncation is t'^/{U/uJo + 2g'^/ujQ). 
Eq. Hll|) gives for Go the expression 



Gn = 



z + 2Mg^/ujo-Ho^i[z- 



Z-^ J 

j 



6|6j + 2(M - l)5V^o - U]2'Hl, 



(12) 



_ffo,i[---]2 ^Hqi is a matrix in a Hilbert space without phonons according to the definition of Pq. Hq^ 
creates phonons as well as a PEDS. This will be picked up by the diagonal matrix [...]2^^. Finally, 
Hq I annihilates the phonons and the PEDS. Consequently, the entire expression is either diagonal or 
has off-diagonal elements with nearest-ncigbor pseudospin exchanges. Such a matrix can be expressed 
by a (anisotropic) pseudospin- 1/2 Hamiltonian. A detailed calculation gives an anisotropic Heisenberg 
Hamiltonian 

Heff = -2MgV^o + Ho,i[z - ^oJ2 + ^(^'^ - l)^'^ - U]2^hI, 

j 

= -2M5V^o + E ['^niS^S?, - i) + auiSfSp + 3^3^,) 
jJ' 

with the pseudospin- 1/2 operators S*^ , S'*' , S"^ and z-dependent coupling coefficients (cf. Appendix A) 

-r r ~ ^""~\ -2gV^i) 



(13) 



.,2 6-'^"/"° ,^U-z-2{M-l)g^/uoQ ^ 2^ 2^ 



Wo 



4 



an = 2< 7 ( /wg). 

Wo Wo 



7* is the incomplete Gamma function 



7 (a,?;) = 2^ 



to! a + TO 

m>0 

Introducing the system energy E ^ z + 2Mg'^/uJo we can write (z — Heff)^^ = {E — Hh)^^, where 
Hjj is the Heisenberg Hamiltonian: 

Hh^Y. [^niS-'S;, - i) + aniSfSf, + SfSy,)\ . (U) 
jJ' 

Now the couphng coefhcients depend on the number of lattice sites M only through E: 

^VA^ r/ + 2gVwo-g 2 2 

an = 2t 7 ( , -2.9 /wo) 

Wo Wo 

OTi =2i 7 ( ,2.9 /wo). (15) 

Wo Wo 

4.1 Crossover from Weak to Strong Electron-Phonon Coupling 

The energy-dependent coefficients in i?e// simplify substantially for weak and strong coupling. Relevant 
are low energies E which represent the poles of the resolvent in Eq. To avoid the poles of the 

incomplete Gamma function, the following discussion is restricted to energies 

E<U + 2g^/LUo. (16) 

It will be shown that in this case there exist poles of the projected Green's function with Ej < 0. The 
restriction (|16|) implies that the approximated diagonal Green's function G2 in Eq. (|ll|l has only negative 
matrix elements: 

G2 = (i? - 2gV^o -U-Y. Ni)^' < 0. 

j 

Thus H01G2HQ1 is a negative matrix and the projected Green's function 

{E - Ho^iG2Hq ^)^^ 

has all poles Ej on the negative real axis. The relevant parameter in our truncated continued fraction 
is U/oJo + 2g^/w§ ^ 1. This allows a free tuning of the electron-phonon coupling g, as long as U is 
sufficiently large. 

The weak-coupling limit corresponds with the Hubbard model. For the latter (i.e. for g = 0) it is 
known that a 1/J7 expansion at half filling gives in leading order an isotropic Heisenberg model with 
coupling coefhcients [SJ 

2^2 

atT = an = -JJ-- 

A similar result was obtained for the Holstein model ^ . It should be noticed that the energy dependence 
disappears in this expansion. For nonzero but small g {g/u>o 1) the coupling coefficients of Eq. (|15|l 
have the asymptotic behavior 

2^2 

an-an^252/^^_^^_^. 

In the opposite regime, where the electron-phonon coupling is strong (i.e. g/ojo 3> 1), the incomplete 
Gamma function is approximated by 



Wo ' ' ^ ' u/ u + 2gyLUo~E ^^m\' U + 2g^/ujo-E- 
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Thus the couphng coefficients of the pseudospin-1/2 Hamiltonian are 

In this limit the effective Hamiltonian is diagonal and accidentally degenerate. 

The crossover regime, within the restriction of Eq. is shown in Fig. 1. It indicates that at weak 

electron-phonon coupling there is a strong isotropic pseudospin-pseudospin coupling, whereas a strong 
electron-phonon coupling implies a strong pseudospin-pseudospin coupling only for the component 
but a weak one for and . Thus the tuning of the electron-phonon interaction is given by a crossover 
from an isotropic Heisenberg to an Ising model. This may be accompanied by a sequence of crossovers 
and/or phase transitions. 

The pole from the groundstate of the projected resolvent {E — Hh)~^ is easily evaluated in the 
asymptotic regimes. For the weak-coupling regime it is 

Eh = ffV'^o + U/2 ~ V(5V^o + U /2f - 2t^XH 

and for the strong-coupling regime 

Ei = g^/LOo + U/2 - Vi9y^o + U/2)^^2tni, 

where Xh (A/) is the lowest eigenvalue of the isotropic Heisenberg (Ising) Hamiltonian with unit coupling, 
respectively. 

5 Conclusions 

Starting from a Hamiltonian with short-range Coulomb and E ® (} Jahn- Teller interaction, a system of 
spinless fermions was studied at half filling. An effective Hamiltonian H^f / was derived under the assump- 
tion that the kinetic energy (i.e. the hopping term) is always dominated by the local interaction energy. 
In the absence of the electron-phonon interaction this leads to the well-known isotropic pseudospin- 
1/2 Heisenberg Hamiltonian for -ffe//- A weak electron-phonon interaction suppresses the pseudospin- 
1/2 interaction of the effective Hamiltonian, and an increasing electron-phonon interaction develops an 
anisotropy, where the pseudospin-pseudospin interaction in the xy plane decreases like exp(— 4(72/a;Q) 
with the electron-phonon coupling constant g and the pseudospin-pseudospin interaction in the z direc- 
tion decreases like 1/g^- 
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Appendix A 



^2^2 V- (0,T |m,0,0)(m,0,0|0,t)(0,i |m^O,0)(m^O,0|0,i) 



(gV^o)'"+'"' 



rn m'>0 "^'"^'^ ^ ~ [^o(™ + "^') ~ {N - 2)g'^ / UJq + t/] ' 



The double sum is reduced to a single sum 

ml z — LOnm 

m>0 

and 



ml z — Lo^m + [N — l^g"^ ju^ — U 
(0, T |m, 0, 0) (m, 0, 0|0, i)(0, j |to', 0, 0)(m', 0, 0|0, T) 



z~[LUo{m + m')-iN ~2)gyu;o + U] 



y ^ 



m!m'! z - [wo(to + m') - (A^ - 2)g'^/ujo + U] 



The double sum is again reduced to a single sum 

{-2gyu;or 



2t2e-29V-o y i 



m! z - cjom + (iV - 2)g^/ujQ - U ' 



These expressions are related to the incomplete Gamma function [7] 

1 i^yr 



such that 



^ — ' ml a + m 

m>0 



7^0 T/ - ^ - 2(M - l)gVc.o 

a^^ = It 7 ( ,-2g /uj„) 

LUq LUq 

e^^9^ t/ - z - 2(M - l).gV^o) ^ , 
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Figure 1: Coupling coefficients of the effective pseudospin-1/2 Hamiltonian of Eq. H14|l as a function 
y — Off (upper set of curves) and (lower set of curves) in units of /uuo for {U ~ E)/loq 

1.0, 1.2, 1.5 (from top to bottom in each set of curves). 
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